Confined nanorods: jamming due to helical buckling 
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We investigate a longitudinally loaded elastic nanorod inside a cylindrical channel and show 
within the context of classical elasticity theory that the Euler buckling instability leads to a helical 
postbuckling form of the rod within the channel. The local pitch of the confined helix changes along 
the channel and so does the longitudinal force transmitted along the rod, diminishing away from 
the loaded end. This creates a possibility of jamming of the nanorod within the channel. 
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I. INTRODUCTION 

Elastic instabilities of nanoscale rods and in particu- 
lar the Euler (buckling) instability have long been rec- 
ognized as being essential for structural as well as func- 
tional aspects of nano- and bio-systems. Some time ago 
[l] Manning proposed that DNA collapse in polyvalent 
salt solutions proceeds through a type of Euler buckling 
driven by diminished electrostatic repulsive interactions 
between charged phosphates along the DNA backbone. 
Though the details of his picture appear to be specula- 
tive, the general framework has been completely vindi- 
cated by later work Buckling instability appears to 
play also an important part in the conformation of linker 
DNA within the chromatin fiber [4] , in the case of a grow- 
ing microtubule pushing against a microfabricated rigid 
barrier [3] and in the stability of AFM-tips, nanotubes 
[5|] and nanorods [6]. 

The confinement of nanoscale rods to micro- and 
nano-channels presents yet another set of experimen- 
tal and theoretical problems [7] centered around the 
technological challenges of manufacturing nanopores and 
nanochannels for investigating and manipulating DNA. 
Here the main theoretical thrust is in the direction of 
understanding the various contributions to the confine- 
ment free energy that depend on the size of the con- 
fining space and the intrinsic properties of the confined 
nanorods. The behavior of confined semiflexible poly- 
mers is particularly important in this context and has 
been recently analyzed in detail [sl, loj . In many respects 
this type of problems are mainly centered upon trans- 
verse confinement, whereas elastic instabilities described 
before are due to longitudinal confinement of elastic rods. 

Motivated by these phenomena we will try in what fol- 
lows to combine the two aspects of polymer confinement 
described above, considering an elastic buckling insta- 
bility within a confining cylindrical channel. The con- 
finement of the elastic rod is thus twofold: longitudinal, 
leading to buckling, and transverse, leading straightfor- 
wardly to a helical postbuckling form. We will show that 
transverse confinement with frictional walls furthermore 
leads to a decay of the longitudinal force within the rod 



along the long axis of the confining channel. This will not 
be too difficult to rationalize since we know from other 
systems that longitudinal stresses can be taken up effec- 
tively by friction at confining walls. For a silo filled with 
granular matter and compressed on one side by a force 
Fo, Janssen's equation [10] gives an exponential decrease 
of the compression force with the distance z from the 
compressed side, 

F{z) = Foexp(-AA:/^z7^/5), 

where A = Orrl (^zz is the ratio of horizontal and vertical 
stresses, V /S \^ the perimeter to cross-section area ratio, 
kfr is the static friction coefficient and we have neglected 
gravity for the purpose. The same holds also for a regular 
elastic solid, i.e., a rigidly confined rod, if one neglects 
lateral deformation gradients in the rod [11,]. There we 
have 

F{z) = Fo ex.-^[-crkfrzV/S{l - a)], 

where <j is Poisson's ratio. On this simple basis one ex- 
pects also a conceptually similar effect, i.e., a decrease 
of compression force with the distance along an elastic 
rod confined to a cylindrical channel in the post-buckling 
regime. This will allow us to hypothesize on the existence 
of a jamming regime for sufficiently long confined elastic 
rods. 

The plan of this paper is as follows: we shall first 
present the classical elastic model for an Eulerian rod 
and derive the scaling relations valid for the onset of the 
buckling instability. We shall solve an approximate form 
of the elastic equations within a cylindrically confined 
channel and show that the longitudinal force, at least 
within the considered approximations, decays inversely 
proportional to the length of the confining channel when 
friction at the cylindrical wall is taken into account. To 
obtain exact solutions of the elastic equations we will 
proceed numerically and derive the minimal shape of the 
confined rod and the magnitude of the transmitted lon- 
gitudinal force along the axis of the confining channel. 
At the end we will discuss what could be the conditions 
that would lead to a jamming of the elastic rod within 
the cylindrical confining channel. 
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II. ELASTIC MODEL 

We model the rod as a thin Eulerian elastic filament to 
which the standard continuum theory of elasticity can be 
applied. The confining channel is modeled as a straight 
cylindrical tube with rigid walls. The ratio R/L of tube 
radius R and undeformed filament length L is the rele- 
vant geometric parameter of the system. 

Let us briefly review the equilibrium equations for thin 
elastic rods fl2|. The force balance reads 



dF 

dT 



K = 0, 



(1) 



where d/ is the length element on the rod (/ is the nat- 
ural parameter — the length along the fllament), F is the 
elastic force exerted on the leading surface of the rod 
element (in the sense of increasing the force on the 
opposite (trailing) side is — F, of course), and K is the 
external force per unit length. The torque balance reads 



dM 
"dT 



-t X F = 0, 



(2) 



where M is the elastic torque exerted about the center of 
the leading surface and t is the unit tangent of the rod. 
The elastic torque is related to the deformation of the 
fllament via 



M Grt^EIt X t, 



(3) 



where E is the Young modulus and we assumed that the 
rod is circular, i.e., both geometric moments of inertia 
of the cross section are equal, / = jdSx^ = fdSy^. 
The flrst term on the right describes torsion which we 
assume to be absent. One can show that in case of the 
isotropic moment of inertia, torsion is absent everywhere 
in the rod when there is no twisting torque applied [12]. 
In what follows, the twisting (torsional) torque will be 
absent, i.e., the torque in the rod will have a vanishing 
tangential component. 



M Elt X t. 
From Eqs. ([2]) and (|4]) one then gets 

Elt xt + txF = 



(4) 



(5) 



as the fundamental equation describing the shape of the 
bent fllament. The solutions of this equation with various 
boundary conditions are discussed in standard references 
on the theory of elasticity [l2[. 



A. Buckling in confined geometry 

It is well known that when a straight rod is compressed by 
an axial force, above threshold it undergoes a buckling 
instability — the so-called Euler buckling — and becomes 
bent. For a rod with hinged ends or ends that are lat- 
erally free (both situations are identical in this case), for 
example, the threshold force is Fcrit = Eln'^/L'^ jl2| . 



What happens when the (thin) rod is conflned within a 
rigid cylinder? Once the compressing force is sufficiently 
large so that the bent rod touches the wall of the cylinder 
(let it be bent in the diameter plane) and is further in- 
creased, there are two scenarios one can think of: the rod 
can bend back-and-forth staying in the diameter plane, or 
it can become a helix touching the cylinder wall (spring- 
like conflguration) . In oil-drilling community it has been 
recognized for a lon g ti me that it is the latter case that 
eventually happens [l3|, [l^ , [TBI . (In fact, this is the 
only area of research to have studied helical buckling.) 
It turns out, as we will show in Section that there in- 
deed exists a secondary threshold above which the planar 
conflguration is unstable with respect to the helix. 



B. Scaling 

By introducing dimensionless quantites denoted with 
tilde. 



EI - EI - 

l = Rl, F=— F, K = — K, 



771 T 

M = —M, (6) 



Eqs. ([I])-® attain a universal form. Thus, the solutions 
depend only on the aspect ratio L/R. Upon rescahng the 
system size and preserving the aspect ratio, the solution 
is unchanged if also the forces and torques of the bound- 
ary condition are rescaled correspondingly, i.e., according 
to Eq. (j6|). Similarly, taking a rod with a different EI 
and rescaling the forces and torques does not affect the 
solution. 



III. ANALYTIC DESCRIPTION OF THE HELIX 

A. Simple model 

Let us write down equilibrium equations for a helically 
deformed fllament constrained to and lying on a cylindri- 
cal surface of radius R. The deformation is sustained by 
an external compressing force F^ < parallel to the long 
axis of the cylinder. We disregard the transition regions 
in the vicinity of both ends and focus on the central part 
of the long helix, i.e., we assume that its pitch is constant. 
Friction is absent at this stage. 

In cylindrical coordinates (r, ^, z) the helix is given by 
r = R and the linear function ^(/). The tangent is thus 
given by 



(7) 



Further we put F = F^Gr + i^^e^ + i^^e^, where the flrst 
two components are unknown and must be determined. 
The force of the cylinder per unit length is K = — i^e^. 
Eq. ([5]) then reads 
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1 - {R^Y{EIR^^ - F^) + R^F, 



l-{R<j>)'^Fre^-R^Fre, 



0. 



(8) 



Eq. ([8]) requires 

while from Eq. ([T]) one gets 

F - ^ 

and hence for the radial component of Eq. 



(9) 



(10) 



1 - {R(j))^ (eIR^"^ + i^) + R^'^F, = 0. (11) 



Now comes a crucial point. The equilibrium condition 
(pTj) involves both the helical deformation (j) and the ex- 
ternal force K, so additional input is needed to determine 
one or the other. This comes from the boundary condi- 
tion, which is however inaccessible under the assumption 
of fixed pitch. 

Let us digress a little and illustrate in physical terms 
why the force of the cylinder, K, cannot be specified until 
the boundary condition is known. Besides the force F, 
the boundary condition involves also the torque exerted 
on free ends (which must be normal to the tangent of 
the rod as we are not considering torsion). By using 
a pair of wrenches and applying torque to the ends, one 
changes the force exerted on the cylinder wall. Using just 
the right torque, for example, one can make K vanish. 
Increasing the torque further, the rod detaches from the 
cylinder and forms a free helix with the radius smaller 
than R [I^. According to Eqs. (g]) and Q, the torque 
in the filament is 



m = EIR(lf 



R^Gz 



(12) 



In fact, with this torque applied to the ends (together 
with and F^), the helical pitch is constant everywhere. 
By changing the torque, one changes the pitch, Eq. (p!2|) . 
and therewith the cylinder force K, Eq. (pT]) . Note that 
there is no tangential force in the free (detached from the 
cyhnder) helix, Eq. (fTQ|) . 

When there is no torque applied to the ends, i.e., when 
compressing the rod inside the cylinder by a pair of pis- 
tons, the part of the rod close to the end detaches from 
the wall while the tip (which is straight due to zero 
torque) pushes against the wall with a discrete force (see 
numeric solutions. Fig. [T]). The torque due to this force 
increases as we move away from the tip and eventually 
becomes sufficiently large for the helical deformation. 

Let us return to the case of a constant helical pitch. 
One can invoke an energy argument to independently es- 
timate the dependence of ^ on the compressing force F^. 




FIG. 1: (color online) A long-axis view of helical filaments 
with lengths (a) L/R = 20 and (b) L/R = 640, buckled by 
a force 1^^(0)1 = 0.25EI /R^ , which is approximately (a) 10- 
times and (b) 10,000-times the critical force for Euler buckling 
Fcrit = F/tt^/L^. There is no torque applied to the ends. 



The elastic free energy of the compressed helix without 
the contribution of the torsion is 



2 io 



dl t 



L,F, 



(13) 



where L is the length of the filament and is the 
length of its projection to the z axis, i.e., the length 
of the coil. Recall that F^ < 0. With Eq. dZj) and 

dz/dl - ./i - ^P^A2 _ 



1 — {R4>)'^ = Lz/L, we have 



LJl-{R4>)^\F, 



and after minimizing, in the lowest order of R4>, 

:o \FJ 



2EI 



(14) 



(15) 



Hence, from Eq. (pTj) one gets a direct relation between 
the load and the wall forces. To the lowest order the force 
of the cylinder is 



K 



R 
2Er 



F\ 



(16) 



So far the friction has been absent. Let us now con- 
sider it as a perturbation and introduce a small friction 
coefficient kfr- One is aware, of course, that for finite 
friction the constant pitch is not a solution, neither can 
one use the energy argument. We are interested in how 
the compressing force IF^I is reduced by the friction go- 
ing from one end of the helix to the other. Let us assume 
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that the hehx is compressed by the external force F^o at 
/ = and that the friction is parahel to opposing the 
external force. Assuming the static friction is maximum 
everywhere, it follows from the force balance ([1]) that 



dl 



kfrK. 



Combining Eqs. (p!6|) and (p!8|) we finally get 

kfrR 



1 



2EI 



The virtual point where would diverge lies at 

le = -2EI/kfrR\F,o\, 



(17) 



(18) 



(19) 



which thus dictates the extent of the reduction of the 
compression force. The smaller |/c|, the strongest is the 
decay of the force. Increasing kfr^ \Fzo\^ oi" R reduces \lc\ 
and moves the singular point closer to the beginning of 
the helix. Eq. ([18]) can be put into a universal form by 
introducing the force unit Fq = 2EI /kfrRL, 



Fo 



Fq 
\F.o\ 



(20) 



Note in Eq. (|6]) that if the aspect ratio is constant, the 
force scales as (or equivalently, l/L'^). This, of 

course, always holds and is approximation-independent. 
In addition, the result ([20]) of our simple model gives 
a 1/RL scaling of the force when the aspect ratio is 
changed, which, however, is only an approximation. 



B. Exact analysis 



For completness, let us derive exact equilibrium equa- 
tions for the helically buckled filament in the presence 
of friction. We do not aim to solve them, yet they will 
be helpful giving us insight into the solution. As before, 
the filament is assumed to lie on the cylinder everywhere, 
r = R, but now all variables are /-dependent, including 
the compressing force Fz. 

Eq. ([5]) involves only two components as it is orthogo- 
nal to t. Therefore we project it to and the direction 
t X Gp. Writing = u;{l)^ the former is 



FIR 



(l-R^uo^) 



3/2 



UJ^ (1 



RWf 



Lj {1- R^uj^) - R^ujuj^ -i^^Vi 



r 



FzRu = 



(21) 



and the latter 



0. 



(22) 



Again we choose from an infinite number of solutions, 
which exist when the static friction is involved, the sim- 
plest and most symmetric one, in which the friction is 
everywhere at its maximum and parallel to z. With 
K = K{l){—er + kfr^z) the three components of the 
force balance ([1]) are 



LoFch — K] e 



(uFr + F^^ ■ 



0, 


(23) 


0, 


(24) 


0. 


(25) 



Eqs. (|2T ]) - (|25]) represent a closed set of five ordinary dif- 
ferential equations for the five variables a;(/), K{1)^ and 
F(/). One can verify that in case a; = 0, Eqs. ([9|)- (pT]) are 
recovered. 

One can expand the system ([2T]) - ([25|) for small Ruj. 
For kfr 7^ 0, terms of third order must be included to 
get a solution (one can check that the second order gives 



only the second order part of the stationary solution ([9])- 
(pT]) for /c/r = 0, and no solution for kfr 7^ 0). Inserting 
Eq. ([22|) into Eq. ([24|) and integrating, one gets a partial 
result F^ = A — EIuo^, where A is a constant. 



IV. NUMERICAL APPROACH 

We will solve the fundamental equations for the actual 
deformation of the compressed and confined filament nu- 
merically for it is only in this way that one can obtain a 
complete solution of the problem. The filament is natu- 
rally allowed to detach from the wall, which essentially 
takes place near the two ends of the filament, and the 
helix within the cylindrical pore can form spontaneously 
without having been put in "by hand". Let us stress 
again that the solution, including in particular the pitch 
of the bulk helix and the force exerted on the cylinder 
wall, is only determined when the ends of the filament 
are taken into account, which furthermore inevitably in- 
volve detachment of the filament from the wall unless 
the boundary condition is very special, as explained in 
Sec. lmAl Eq. (p^. 

For the purpose of numeric modelling that involves dis- 
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cretization, it is generally better, if only possible, to start 
with a discrete analogue of the continuum system and 
write down algebraic equations, rather than discretizing 
the differential equations themselves, derived for the con- 
tinuum limit. In this spirit, the elastic filament will be 
represented as a set of straight and stiff elements (links) 
of fixed length Iq with forces and torques acting between 
them. The element i is described by its center of mass 
and a unit (tangent) vector giving its orientation. 
Dynamical evolution of the filament shape is obtained 
by Newton's laws for translation and rotation of each 
element. The joints between the links are rather cum- 
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FIG. 2: The external load dependence of the projection of fil- 
ament length onto the axis of the confining cylindrical pore. 
The weak force dependence of length in the straight configu- 
ration is due to the small artificial compressibility introduced 
by the penalty potential (|26)) (i.e., finite value of the spring 
constant k). In the limit of zero compressibility this part of 
the functional dependence would be a straight horizontal line. 



R) 



R 
R 



(27) 



Again the coefficient should be sufficiently large so 



that 



-i? <C lo- The forces on the ends of the elements 



dV 



dr, 



± ' 



(28) 



It makes little difference (i.e., no difference in the limit 
lo 0) whether the force of the wall actually acts on 
both ends or just at the center of the element. 

Torques about the center of element i are of two kinds. 
One comes from elastic couples exerted by the two neigh- 
bouring elements and the other from the forces on both 
ends of element i: 

M, = Ct, X (t,_i + t,+i) + ^ (t, X F+ - 1, X Fr) , 

(29) 

where C is the bending stiffness which we shall connect 
to continuum parameters. In the continuum picture, the 
first part of the torque translates to the first term of 
Eq. ([5]), and the second part to the second term of Eq. ([5]). 

To get the evolution of the shape of the filament, we 
coveniently assume overdamped dynamics: 



dr^ _ I _ 

dt ^ ^ ^ ' 
dt^ d(p^ 



dt 



dt 



M. X t., 



(30) 
(31) 



where (3 is an arbitrary damping coefficient defining the 
time scale. 

One can verify that the set of discrete equations ([29])- 
(|3T]) agrees with properly discretized continuum equa- 
tions with the following connection between the param- 
eters: 



bersome to model: they act as constraints for the ends 
of the two elements which must meet in a single point, 
and thus exert constraint forces that can be determined 
only implicitly, such that the elements satisfy the con- 
straints. To avoid this, we relax (soften) the constraints 
and introduce a quadratic (bond) penalty potential 

V{vt,vT^,) = ^\vT^,-vt\\ (26) 

where r+ = + (/o/2)ti and r^:^^ = r^+i - (/o/2)ti+i 
are the two ends meeting at the joint. The coeffi- 
cient k should be sufficiently large so that the gap be- 
tween the ends is small compared to the element length, 
|r^-L — r^l <C Iq. The side wall of the cylinder is mod- 
elled in a similar manner, introducing a wall penalty po- 



C = (32) 

while (5 = lo/3ph^ where /3ph is the physical (continuum) 
damping coefficient. The latter is only important for 
properly scaling the time when changing the element 
length lo in the numerical model. Similarly, k = kph/lo 
to ensure that the spring constant kph of the rod (which 
is infinite in the continuum description ©-([I]) of a thin 
rod) is unaffected by changes of Iq. 

To avoid giant friction forces in certain discrete points 
(e.g., a few points arround the detachment region which 
exert a large force on the confining wall), the friction 
force is properly capped, so that the Coulomb's friction 
law is violated in these few points but met everywhere 
else. 
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V. RESULTS 

First let us take a look at numeric solutions of heli- 
cally buckled filaments. Fig. [T] shows two examples of 
filaments of different lengths buckled by the same axial 
force. We use the most natural boundary condition: the 
upper (compressed) end is only pushed along z (the cylin- 
der axis), while the lower end is not allowed to move in 
z direction. Both ends are free to move laterally until 
they hit the confining wall. There is no torque applied 
to them. To build up the torque exerted on the filament 
cross section in the helical state, however, the ends must 
leave the helical state as anticipated in Sec. IIII A[ 




FIG. 3: (color online) Planar configuration of the confined 
filament (a) is unstable (b) and is transformed into the heli- 
cal (c) upon perturbation. The helix can contain met ast able 
irregularities as seen in the equilibrium state (c). 



To numerically study the solution path from the 
straight to the helical configuration, we introduce a tiny 
random perturbation to the straight compressed fila- 
ment. A numerically obtained state diagram is depicted 
in Fig. [21 where a suitable parameter distinguishing be- 
tween the configurations is simply the projection of the 
filament length onto the z axis. There exist two threshold 
forces beyond which the straight and the Euler-buckled 
states are unstable, respectively. The helical state is a 
consequence of the confinement and sets in after the fil- 
ament has touched the wall. Above the threshold, the 



planar state remains a solution but is unstable with re- 
spect to the helix, which is demonstrated in Fig.[3l There 
one also sees a defect in the helix. Such defects are a 
general feature subject to the initial condition (pertur- 
bation) and are metastable, i.e., they are stable while 
the compressing force is not reduced. 
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FIG. 4: (color online) Compression force profile |i^2(OI ^ 
filament with length to radius ratio 320. The loading force is 
\Fz(0)\ = 0.25^//i?^ which is approximately 2,600-times the 
critical force for Euler buckling Fcrit — Eli{^ jl? . Numeric 
data (black) is excellently fit (red) with the function /(/) = 
+ /)^ (the curves are perfectly overlapping on the scale 
of the figure): (a) A = 36.3, /o = 285.5, c = -0.882, (b) 
A = 21.4, /o = 131.8, c = -0.894, (c) A = 11.95, /o = 44.9, 
c = -0.903. 



We are furthermore interested in the force that has to 
be applied to the lower end to keep it fixed, i.e., to sustain 
the confined helix. As we are to show, this force can be 
much less than the load applied to the upper end due to 
the action of friction and in fact approaches zero for very 
long confining channels. This opens up the possibility 
of jamming, i.e. of a stable static helical configuration 
of the rod that has been jammed against the confining 
walls via the surface friction and is sustained with only 
a tiny (zero in the limit of very long channel) opposing 
force at the other end of the channel. 

Fig. m shows the profiles of the longitudinal force IF^I 
as we move along the filament, and Fig. [5] shows the pro- 
file of the helical deformation (\). The force profiles are 
perfectly fit by a power law, yet with an exponent close 
to -0.9 instead of -1 as suggested by Eq. (p!8|) . The va- 
lidity of the power law is remarkable though it deviates 
from the approximate universal scaling of Eq. ([2Q|) . In 
other words, the solution is practically indistinguishable 
from a power-law, in spite of the evidence that it con- 
siderably departs from our simple model. The reason for 
this remarkable validity of the power-law, yet with an 
exponent different than -1, is not clear. Fig. [6] shows the 
force transmitted to the lower end as the function of the 
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FIG. 5: Profile of d(/)/dZ, the inverse of tiie iielical pitcii qo = 
27r/(d(/)/d/), for tiie filament of Fig. |4] (c). The boundaries 
where the pitch is not defined are not displayed. As expected, 
(j) gets smaller (the pitch gets larger) as the longitudinal force 
is reduced. Its relative decrease is however smaller than that 
of the force, as hinted by the quadratic dependence Fz oc (p^ in 
Eq. ([15)) of the simple model (which assumes constant pitch). 



simple model. On the other hand, increasing the loading 
force beyond a threshold (that scales as results in a 

catastrophic event — a U bending of the filament, leading 
to the escape out of the confining pore. Fig. [HI Hence, it 
is only for long and thin rods that the transmitted force 
can be strongly reduced. 
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length L of the filament, at a fixed load on the upper 
end. We see that for long filaments the transmitted force 
is strongly reduced and asymptotically approaches zero. 
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FIG. 6: Compression force, \Fz{L)\, transmitted through 
filaments of lengths L; the external compressing force is 
IF^O)! = 0.25EI/R'^; kfr = 0.8. The dotted line serves as an 
eye guide. 

The physically interesting interval of loading force 
strengths depends on the aspect ratio of the filament. In 
Fig. [31 for example, the loading force was selected such 
that the transmitted force was significantly reduced. The 
reduction of the transmitted force is weaker if the load- 
ing force is smaller. Fig. [71 in accord with Eq. ([T9|) of the 



FIG. 7: (color online) Compression force profile |i^2(/)| (black) 
in the filament with length to radius ratio 320 for a smaller 
loading force \Fz(0)\ = 0.0125^//i?^ and kfr = 0.8, fit (red) 
with the function /(/) = A{lo + /)^; A = 15.5, lo = 2010, 
c = —0.942. The reduction of the force is much smaller than 
for larger loads in Fig. [4] on account of the singular point Ic in 
Eq. (|19p moving further away. The exponent of the power law, 
however, remains nearly unchanged. The boundary regions 
get wider as the load is reduced. 



VI. DISCUSSION 

In this work we have analysed the behavior of an elas- 
tic filament that is confined in the longitudinal as well 
as transverse direction. The longitudinal confinement 
leads to the well known Euler buckling instability that 
depends on the external loading and the length of the fil- 
ament. We have shown that the buckled configuration on 
touching the confining walls of a cylindrical channel then 
evolves through a planar deformed configuration towards 
a helical state via an additional instability. Considering 
the effects of the friction on the walls of the cylindrical 
enclosure we have been able to demonstrate, that the 
longitudinal force transmitted through the filament de- 
cays along its length. This phenomenon is not unrelated 
to the decaying longitudinal stresses within a cylindrical 
granular column anchored by the wall friction, or even a 
regular elastic solid with a finite Poisson ratio enclosed 
within a rigid hollow cylinder and again anchored by wall 
friction. 

The important difference between the examples of 
granular and regular solids confined within cylindrical 



8 




FIG. 8: (color online) If the load is too large, the filament 
makes a U-bend and escapes out of the cylinder as shown on 
the subsequent figures instead of forming the helix. L/R — 
320, |Fz(0)| = 0.5EI/R'^, only the upper half of the cylinder 
is shown. 



walls and the present case of thin elastic filament is the 
nature of the decay of the longitudinal force along the 
cylindrical enclosure. In the first two cases the decay is 
exponential and leads to a natural length scale for the 
problem. This means that the jamming of the granular 
or classical elastic bodies depends only on intrinsic pa- 
rameters, describing the stress distributions within the 
body and the magnitude of the friction forces between 
the body and the cylindrical enclosure. In the case of 
the confined filament the longitudinal force transmitted 
along the cylindrical enclosure decays algebraically. In 
this case the onset of jamming is determined by the in- 
trinsic elastic parameters of the filament as well as the 
longitudinal force counteracting the loading from the op- 



posite side of the enclosure. This counter-force could be 
in principle very small and could even result from ther- 
modynamic fluctuations on the other end of the enclosure 
if the cylindrical pore is small enough. Nevertheless the 
fact that the onset of jamming in the case of a confined 
elastic filament is scale free and thus depends on external 
constraints separates it fundamentally from the standard 
jamming in granular materials [18]. 

What would be the systems that could exhibit this 
type of jamming scenario? We have found out that the 
criterion for the onset of helical postinstability shape is 
geometrical: according to the curve in Fig. [6l the as- 
pect ratio of polymer length (or persistence length, if it is 
smaller) and the radius of the pore should be of the order 
of 100 or more for a significant reduction of the compres- 
sion force that could in principle lead to jamming. For 
carbon nanotubes with a typical length of 100 jam and 
Young's modulus of 10^^ Pa, the relevant pore radius is 
1 /im or less, and the loading force is in the nN range. All 
reasonable values. On the other hand, DNA appears to 
be too flexible, i.e., its persistence length of 50 nm is too 
short as it should be confined to a sub nanometer pore, 
violating the structural integrity if DNA (the diameter 
of DNA is approximately 2 nm) . It thus appears that at 
least in principle one could observe the helical instabil- 
ity coupled to a jamming transition within the context 
of confined microtubules. We hope that our theoretical 
work will provide enough motivation for experimentalists 
to search for this interesting phenomenon. 
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